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TECHNICAL MEMORANDUM 


A FIELD STUDY OF SOLID ROCKET EXHAUST IMPACTS 
ON THE NEAR-FIELD ENVIRONMENT 

L INTRODUCTION 

The launch facilities for the Space Transportation System (STS), the space shuttle, 
are unique in the sense that copious quantities of water are sprayed into the base region 
below the vehicle and into the flame ducts during each launch for the purpose of 
dampening the acoustic and initial overpressure waves generated by the vehicle. At the 
Kennedy Space Center (KSC) launch facilities, this amounts to approximately 900,000 
liters from when the flow begins (a few seconds prior to main engine ignition) to 10 sec 
after lift-off. If the shuttle launch facilities at Vandenberg Air Force Base are acti- 
vated, the flow is expected to be roughly twice as great. It has been shown (1) that this 
amount of water exceeds the amount that can be vaporized by the exhaust heat. Some of 
the excess is atomized, mixed with the exhaust, and results in a deposition of hydro- 
chloric acid and solid material (mostly aluminum oxide from the solid rocket motors). 

The deposition is quite heavy near the launch pad. In trace amounts it has been detected 
as far as 22 km downwind on some occasions. Its impact has been characterized by a 
number of careful studies conducted at KSC (2-4). The chemistry and ice nucleation 
properties of the solid fraction of the deposition are discussed in Refs. (5)-(7). 

The properties, location, and behavior of this deposition are of interest primarily 
because of the potential for impacting launch pad operations and the near-pad environ- 
mental quality. The acid content is great enough that it impacts both vegetation (8) and 
animal life near the launch facilities (9). Likewise, it can be highly corrosive to man- 
made structures. It has also been found (10) that gaseous hydrogen chloride is released 
into the atmosphere as the deposition dries. This "revolatilization" process is of poten- 
tial concern in both the areas of human health and corrosion control. 

This report presents results from a series of field studies and analysis which were 
undertaken to help quantify and understand the near-field effects of this deposition. 
Primary emphasis was given to measuring and understanding the effects at KSC launches 
in order to provide a basis for developing reasonable estimates of what to expect from 
launches involving either new vehicles or new launch facilities. Measurement and analy- 
sis techniques suitable to the situation were developed. The work developed out of 
previous studies of the far-field effects of rocket exhausts (1). It was conducted over a 
period of more than 4 years and included field studies at two shuttle launches, 41 D and 
51 A. Additional studies of the exhaust cloud properties from these two launches are 
reported in Refs. (11) and (12). 

An important aspect of this work was the use of the Acoustic Model Facility of 
the Test Laboratory at the Marshall Space Flight Center (MSFC). Here a 6.4 percent 
scale model of the space shuttle is statically fired to study the acoustic environments 
and other phenomena produced during a launch. From May 1982 through April 1984, a 
series of test firings were conducted to examine the initial overpressure (IOP) wave to be 
expected from Solid Rocket Booster (SRB) ignition of the space shuttle at the 
Vandenberg Air Force Base launch site (VLS), the acoustic environment of a VLS shuttle 
launch, and the acoustic environment of the aft cargo carrier in a KSC launch. Selected 
tests from among these were monitored to study the production and properties of acidic 
deposition produced during a shuttle launch. This provided an important opportunity to 



develop and test measurement techniques and to study deposition formation in the VLS 
configuration. 


IL BACKGROUND INFORMATION 
A. Model Description 

In the Acoustic Model Facility the shuttle SRBs are modeled by Tomahawk Missile 
solid rocket motors manufactured by Morton-Thiokol Chemical Corporation. The motor 
contains 175.57 kg of type TP-H-3095 propellant (20.4 percent aluminum) which burns in 
approximately 9 sec. Chemically this propellant is very similar to the propellant in the 
shuttle SRBs. The average rqass flux from the motor is thus 19.5 kg s _ * which is a factor 
of about 4.096 x 10 -3 = 0.064" smaller than the typical mass flux from the shuttle SRBs. 
Actually the mass fluxes of both the model and shuttle motors vary with time during the 
burn cycle. Likewise, propellant temperature and changes in shuttle SRB design also 
make a difference. The shuttle began using higher performance motors with STS-8. 

For the purposes of this study, a value of 5560 kg s~* is used as representative of 
output from each shuttle SRB. This figure represents the average mass flux in the first 
18 to 20 sec following ignition as computed for the preflight analysis of mission STS-13. 
The output is less later in the burn. Before STS-8 the output was about 5400 kg s , a 
value which is still within the ±5.3 percent variation which may occur because of changes 
in the temperature of the propellant. 

The scaling of the model for acoustic and initial overpressure studies is based on 
the ratio of mass fluxes from the solid motors. Linear dimensions are scaled by 0.064 so 
that areas are scaled by 4.096 x 10~ 3 (0.064 squared); thus the model is commonly known 
as the "6.4% Model." As illustrated in Figure 1, the launch mount including the flame 
trenches and launch platform is modeled in plate steel and "Fondu Fyr," a concrete-like 
refractory material. The orbiter and external tank are also modeled to the same scale. 
The orbiter model contains three working engines fueled by gaseous hydrogen and liquid 
oxygen which were originally used to model the Saturn J-2 upper stage engine. The 
scaling factors for the model are summarized in Table 1. 

The most important parameters involved in the formation of acidic deposition in 
the launch process are the mass of HC1 released, the thermal energy released, and the 
volume of cooling water in the cloud/flame trench system. Note from Table 1 that the 
exhaust mass flow rate is scaled by the factor 0.064 2 = 0.0041. Since propellant compo- 
sition in the model is very similar to the shuttle propellant, the HC1 and the thermal 
energy fluxes scale very closely to the same factor. To estimate the relative total mass 
of HC1 and thermal energy, the interaction time between the shuttle SRM plume and the 
on-pad water must be known. From launch photographs, the interaction time is esti- 
mated to be 7 to 10 sec, compared to the 9 sec that the Tomahawk burns. On the model, 
IOP /acoustic suppression water flow rate is also scaled by the same 0.064 2 factor. 
However, for tests in the VLS configuration, the water flow time interval prior to igni- 
tion is much shorter in the model, less than 1 sec, compared to the 15-sec full scale. 

This timing difference is necessary to maintain the properly scaled cross sectional area 
in the duct. As a result, the total volume of cooling water in the model VLS system is 
scaled by a factor much smaller than 0.064 . Water flow into the ducts is illustrated by 
Figures 2 and 3; scale factors are given in Table 2. 
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Figure 1. "6.4 percent model" facility configured for Western 

Test Range testing. 
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Figure 2. Total cumulative water flow per SRB at Kennedy Space Center which may interact with 
the solid motor exhaust. (Flows on the SSME side are excluded.) 
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TIME (SECONDS) 

Figure 3- Total cumulative water flow per SRB planned for the Vandenberg launch site which 
would interact with the solid motor exhaust. (Flows on the SSME side are excluded.) 




TABLE 1. COMPARISON OF MODEL AND FULL-SCALE* PARAMETERS 


Parameter 

Full-Scale 

n 

Full x 0.064 n 

Model 

Shuttle SRM (one) 

Sea Level Thrust (Newtons) 

13.3 x 10 6 

2 

54,700 

47,200 

Mass Flow Rate (kg s~*) 

5560 

2 

22.8 

21.2 

Exit Area (m^) 

11.3 

2 

0.0465 

0.0406 

Expansion Ratio 

7.72 

0 

7.72 

7.36 

Exit Half Angle (deg) 

6 

0 

6 

15 

Exit Mach Number 

2.95 

0 

2.95 

2.81 

Exit Diameter (cm) 

380 

1 

24.3 

22.7 

Supersonic Core Length (m) 

59 

1 

3.77 

3.35 

SSME (one) 

Sea Level Thrust (Newtons) 

1.68 x 10 6 

2 

6870 

6790 

Mass Flow Rate (kg s -1 ) 

472 

2 

1.93 

1.81 

Exit Area (m^) 

4.17 

2 

0.017 

0.0062 

Expansion Ratio 

77 

0 

77 

8.5 

Exit Half Angle (deg) 

5.4 

0 

5.4 

20 

Exit Mach Number 

4.234 

0 

4.234 

3.325 

Exit Diameter (cm) 

230 

1 

14.7 

8.9 

Supersonic Core Length (m) 

54 

1 

3.5 

1.6 


♦SSME at 100 percent power, high performance SRM. Data are reduced to three 


figures since motor performance varies several percent with time and temperature. 


B. Background on Acid Deposition Production 

The data presented in Table 2 illustrate a critical point relative to formation of 
deposition from the STS exhaust; an excess of water enters the pad/flame trench system 
over and above the amount that can be vaporized by the available exhaust heat. The 
excess, which is of order 200 kl (50,000 gal) per SRB at KSC, 500 kl (130,000 gal) per SRB 
at Vandenberg, interacts with the exhaust plume. Much of this water is atomized by the 
mechanical shears and turbulence generated by both the Space Shuttle Main Engine 
(SSME) and SRB exhaust flows. It is then expelled into the near field or mixed into the 
exhaust "ground cloud," which lifts and blows away with the wind. In this interaction, it 
scavenges significant quantities of SRB exhaust products; gaseous hydrogen chloride and 
aluminum oxide particles. The composition and spatial distribution of this material is of 
primary interest in this study. 
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The IOP/acoustic suppression water is injected around the launch mount "holes" at 
the base of the vehicle. In the KSC configuration, water is also injected into the flame 
trenches and spread out onto the surface of the mobile launch platform. In both configu- 
rations, most of the interaction between water and exhaust occurs in the exhaust ducts. 
The excess water is ejected with the exhaust. At Kennedy, both SRBs are directed into 
one duct which empties horizontally toward the north. There are no structures to deflect 
the plume upward, but thermal buoyancy causes lifting after the initial jet decays. The 
KSC SRB exhaust duct configuration is shown in Figure 4. 

At Vandenberg the situation is quite different. First, each SRB is directed into a 
separate duct, one toward the north and the other south. Second, the ducts are covered 
and of a more complex design, as illustrated by Figure 5. The covered portion of the 
north duct runs horizontal for 16 m, then upward at an 8-degree angle for 40 m. For the 
south duct, the horizontal run is 21 m, and the 8-degree run is 35 m. Both ducts open 
onto an uncovered apron area which extends upward at 15 degrees to ground level. 

Beyond the exit of the duct on the north side, the ground level drops gently beginning a 
few tens of meters beyond the opening. For the duct opening to the south, it drops 
quickly into the side of a small ravine, then rises fairly abruptly on the far side to a level 
above the top of the apron. The crest on the far side is less than 300 m from the launch 
mount; its elevation is such that it lies on the projection of the centerline of the 8- 
degree portion of the duct. 


ffl. FIELD TEST RESULTS 
A. Observations of Model Firings 

The spray pattern of the acidic deposition in the immediate vicinity of the launch 
mount is primarily dictated by the exhaust duct and launch mount design. The projection 
angle in the vertical for the material ejected from the north duct was measured on the 
March 25, 1983, model firing. This was a 100 percent of baseline water flow, VLS con- 
figuration test with the vehicle at the zero elevation level (resting on the pad). An 8-m 
vertical pole was placed 31.5 m from the vehicle directly in line with the north duct. 

Test tubes were taped to the pole every half meter, each tilted toward the exhaust at 
about a 45 -degree angle. Figure 6 illustrates the amount of material deposited in the 
tubes as a function of the projection angle above the horizontal relative to the lip of the 
apron. The maximum amount of material was found in the 8- to 10-degree range, indi- 
cating that the covered portion of the duct controls the elevation angle more than the 
15-degree open apron. The distance of the measurement site to the vehicle, 31.5 m, 
corresponds to nearly 500-m full scale (applying a simple 0.064 scale factor). 

In more than half of the 20 test firings that were studied, the primary objective of 
the measurement effort was to determine the deposition pattern on the ground. All of 
these tests were of the VLS configuration with 100 percent of baseline water flow. In 
some cases, the vehicle was elevated above the pad to simulate lift-off. Deposition 
measurements were made by setting out an array of collectors in the target area and 
volumetrically measuring the amount of deposition collected. In the first test, plastic 
petri dishes and 1 -liter beakers were tried as collectors. The petri dishes were blown 
away by the exhaust and so were not usable. Beakers worked fairly well, but some of the 
exhaust material remained on the side walls rather than running down into the beaker. 
Even though the time from rocket firing to analysis was minimized (typically 1 hr) there 
was still concern that a substantial fraction of the deposition might evaporate before the 
beakers were gathered and the contents measured. Therefore, plastic cups with sealable 


8 


SRB TRENCH SSHE TRENCH 



9 












Relative amount of Deposition Collected 


Deposition Collected on Pole 31.5m from Model. 
WTR Baseline Configuration 


3.0 


2.0 


1.0 


* 1 


8 10 

Angle above Lip of Trench 


12 


14 


Figure 6. 


Angular distribution of deposition measured from the horizontal in a 
Western Test Range (Vandenberg) model configuration. 
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lids, 8.1-cm diameter by 6.3-cm deep, filled to within about 1 cm of the top with silicone 
oil, were used. The oil cups were much more labor-intensive and messy, but they pro- 
vided more reliable results. Comparison with the first tests using beakers shows es- 
sentially identical results with both methods. The beakers, however, were used only 
during February and March when the temperature was 5 to 13 °C and evaporation was 
slow. The comparison would probably be less favorable if the beakers were used in 
summer. 

Figures 7 through 10 display four of the measured deposition patterns, two with 
calm or very light winds and two with moderate winds. The outside contour shown in 
these figures represents the limit of the area completely wetted by deposition, or about 
0.1 mm of liquid depth. This contour is easily located in the area to the left-hand side of 
the duct centerline because this area is paved out to 80 m from the vehicle and perma- 
nently marked with a polar coordinate grid. Immediately after a firing, the full wet/ 
partial wet boundary was marked in chalk. However, 3 to 4 m to the right-hand side of 
the duct centerline, the pavement ends and the area is grass covered. A broad ditch, 5-m 
wide by 1.5-m deep, parallels the duct centerline with its center displaced about 12 m to 
the right. In this area, and beyond 80 m in all directions, the full wet/partial wet bound- 
ary must be estimated from the collector cup data with some loss in precision. 

With little or no wind and the model vehicle at the zero elevation level (resting on 
the pad), the fully wetted area extends to about 80 m from pad center. Maximum 
breadth occurs at a range of about 50 m and is of order 35 m. To interpret this as an 
expectation for full scale one must apply a scale factor which, for the sake of simplicity, 
we assume to be 0.064. This simple linear scale factor would be expected to apply for a 
single phase, constant density gas jet. In actuality we are dealing with a two-phase flow 
with substantial temperature differences. Thus drag, gravitation, and buoyancy may be 
expected to play an important role. The 0.064 factor will lead to an over estimate of 
effects in the far field. Scaling the low wind values given above implies an upper limit 
range of about 1250 m and a maximum width of about 550 m at 800 m from the pad. 

With the vehicle raised above the pad the pattern narrowed and lengthened. At the 
maximum vehicle elevation tested, equivalent to 36.6-m (120 ft) full scale, the fully wet 
area extended to about 150 m which corresponds to a 2300-m full scale upper bound, 
measured from the pad center. As a point of comparison, observations at KSC indicate 
that the fully wet area extends 600 to 700 m north of the pad center. From a single 
observation made with the model in the KSC configuration, it appears the 0.064 factor 
over estimates the length of the wet zone by a factor between 1.3 and 2.1. For rough 
"best estimates" a scale factor of (1.7 x 0.064) = 0.11 is probably suitable. 

The data set was also analyzed to obtain an indication of the range from the 
vehicle where the maximum depth of deposition may be expected. The location of this 
point varied from test to test. The lowest value observed occurred on the March 9, 1983, 
test (Fig. 9), which was certainly influenced by the wind. The value of 15 m would scale 
(x 1/0.11) to about 140 m at Vandenberg. The maximum distance of 40 m was observed 
on the September 9, 1983, test in a 4 m s -1 crosswind (Fig. 8). In this case the plume 
centerline was blown to the side of the main array of collectors so the true maximum 
may have been missed. At full scale, 40 m corresponds to about 380 m. The best data 
set for the vehicle at zero elevation (the March 25, 1983, test, illustrated in Fig. 7) and 
one other poor set show the maximum at 20 m (180-m full scale). Two tests with the 
vehicle elevated slightly gave maxima at 30 m (270-m full scale). However, since better 
than half of the water which becomes deposition is already in the trenches at L = 0 for 
both the Kennedy and Vandenberg launch sites, the actual full-scale pattern will tend to 
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MARCH 25, 1983 MODEL TEST 
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Figure 7. Deposition pattern for the March 25, 1983, VLS configuration model test. Depth figures 

are for total deposition, liquid plus solid. Distances indicated are for the model system. 








DEPOSITION DEPTH (MILLIMETERS) 
SEPT. 9, 1983 MODEL TEST 
WTR BASELINE 
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DEPOSITION DEPTH (MILLIMETERS) 
MARCH 9, 1983 MODEL TEST 
WTR BASELINE 
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Figure 9. Deposition pattern for the March 9, 1983, VLS configuration model test. Depth figures 

are for total deposition, liquid plus solid. 



DEPOSITION DEPTH (MILLIMETERS) 
SEPT. 1, 1983 MODEL TEST 
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Figure 10. Deposition pattern for the September 1, 1983, VLS configuration model test. For this 
test, the model was anchored in an elevated position to simulate lift-off conditions. 


look like the low vehicle elevation patterns. Thus at Vandenberg the maximum deposi- 
tion depth may be expected at the 200 to 360 m range for calm or light crosswind launch 
conditions, assuming a 0.11 scale factor. 

The deposition from both the model and actual shuttle launches is a solid/liquid 
mixture; hydrochloric acid, aluminum oxide granules from the solid rocket motors, bits of 
sand, seeds and other debris from the near pad area, and various dissolved trace ele- 
ments. For the September 1 and September 9, 1983, model tests, the solid fraction of 
the deposition was found to be 11 percent in both tests with standard deviations of 4 and 
5 percent, respectively. (The smallest samples were disregarded in obtaining these 
figures because the measurement is not believed to be reliable when the total amount of 
material is less than 0.5 ml. In the very small samples the solid fraction tends to appear 
much larger. Their inclusion in the data would drive the averages to 18 and 13 percent 
with standard deviations of 16 and 10 percent.) By comparison, samples from actual 
shuttle launches show a greater solid fraction. A single sample collected under oil during 
the STS -4 launch was 30 percent solid, 70 percent liquid. During the launch of STS 41 D, 
samples were taken from an array of oil dishes identical to those used for the model 
tests. The average solid fraction for 12 samples was 27 percent with a standard devia- 
tion of 10 percent. 

Apparently the solid fraction is determined by the rate of scavenging and the 
amount of time the drops spend in the cloud. At full scale the drops have much more 
time to collect solid material. At Vandenberg the physical size of the cloud should be 
about the same as at KSC, so one would expect about the same solid fraction in the 
deposition, about 30 percent, unless the aluminum oxide supply is limited relative to the 
amount of water; the amount of water is greater at VLS. 

The data on the acid content of the deposition samples collected from the various 
model tests can be summarized as follows: 

1. Eleven samples collected in beakers, without oil, in four test firings, all on cold 
days. Temperature at test time ranged from 5 to 13 C: mean acidity = 1.26 normal (N), 
std. deviation = 0.20. 

2. Samples collected under silicone oil (Dow Corning 200 Fluid, 500 cs) in 
"Freezettes," 9-cm diameter cups, from a test on a warm day. Model at "zero level" 
elevation, resting on the pad. Fifteen samples from September 9, 1983, test; tempera- 
ture = 28 °C: mean acidity = 1.69 N, std. deviation = 0.27. 

3. Samples collected by same method as 2 above, but from a test with the model 
elevated to the "30-ft" level. Thirteen samples from August 26, 1983, test; temperature 
= 32 °C (see Fig. 11): mean acidity = 1.64 N, std. deviation = 0.39. 

4. Samples collected on warm days by the "milk stand" method. The bottom was 
cut out of a cleaned, polyethylene, gallon-size milk container that was mounted in an 
inverted position over a sample bottle. The milk container formed a large funnel and the 
sample jars could be quickly capped after the firing. The results are summarized in 
Table 3. 

This data set contains considerable variability which is traceable to several 
causes. Some samples taken on or near the 0-degree azimuth show a systematically 
lower acid concentration than those from the 180-degree azimuth. Certainly this may be 
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expected because of the additional water from the SSME duct that is exhausted in this 
direction. It was also observed when examining the deposition drops in the oil cups, that 
some of the drops were clear liquid (plus solids) and some were definitely a light green 
color (plus solids). Both clear and green drops were found together in the same cups. It 
is suspected that the color difference is due to a chemical reaction of the IIC1 with some 
substance in the solid fraction of the deposition which is randomly mixed into some drops 
and not others. The effect is seen in samples collected well away from the influence of 
the SSME duct exhaust and it is not the effect of dew fall since the oil cups were usually 
deployed in the late morning or afternoon, a few hours before the firing. 


TABLE 3. DEPOSITION NORMALITIES FROM "MILK STAND" 

COLLECTORS 

Test Date 

Temp. (°C) 

Location 

Normality 

Aug. 26, 1983 

32 

180°, 

25 

m 

2.07 



180°, 

30 

m 

1.94 



180°, 

40 

m 

2.63 

Sept. 1, 1983 

33 

0°, 

30 

m 

1.40 



0°, 

50 

m 

1.56 



10°, 

30 

m 

0.75 



180°, 

40 

m 

1.63 

Sept. 9, 1983 

28 

0°, 

25 

m 

1.70 



o°, 

35 

m 

1.83 



0°, 

45 

m 

2.19 



180°, 

25 

m 

2.20 



180°, 

35 

m 

2.40 



180°, 

45 

m 

2.45 

Sept. 15, 1983 

27 

O 

o 

00 

rH 

30 

m 

1.84 



180 , 

40 

m 

1.97 



180°, 

55 

m 

2.05 


Samples collected with the "milk stands" tend to show a higher acid concentration 
than samples collected under silicone oil. This is probably the result of water evapora- 
tion from the funnel surface; note that the "milk stands" were only used in the warm 
weather tests. To verify the silicone oil method, 20-ml samples of 2.5 N HC1 were left in 
sample bottles for 10 weeks, one tightly capped and the other under 2.3 cm of oil (open 
part of the time). When titrated, the acid concentrations were equal within the expected 
accuracy of the titration, indicating that significant amounts of HC1 were not absorbed 
into the silicone oil. Acid concentration of the cold weather samples collected in 
beakers is lower than the warm weather samples, indicating a significant temperature 
effect. The beaker samples must have been concentrated somewhat by water evapora- 
tion in the 30 to 60 min required to gather and transfer the samples into closed bottles, 
although at the cold temperatures involved the effect may not have been too great. 
However, the mean concentration, 1.26 N, is still less than the samples collected under 
oil on warm days, 1.64 and 1.69, but the difference is about the same magnitude as the 
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standard deviations in the data sets. Thus the temperature effect cannot be precisely 
quantified because of the scatter in the data. The temperature effect is even more 
noticeable when comparison is made with the "milk stand" samples which averaged 2.12 
N for samples from the 180-degree azimuth. 

Looking at Figure 11 and the data tabulated in Table 3, it is also noted that there 
is a systematic increase in acid concentration with range from the model within each 
data set. If the 11 data points from the 4 cold day tests are plotted together, they 
likewise show an increase in normality with range. In all cases, the scatter in the data 
set is large compared to this effect; but, considering the complete data set, it appears 
that the increase is of the order of 0.02 N/m. 

B. Shuttle Launch Observations 

A small field program designed to verify HC1 revolatilization was conducted at 
KSC in conjunction with two launches in the late summer and fall of 1984. A time- 
resolved, spatially-averaged measurement of the total HC1 concentration was desired to 
provide a general picture of the post-launch work environment on the pad. Hydrogen 
chloride can occur as either gas or as aerosol in combination with water, so simultaneous 
measurements of both forms are required. Both gas and aerosol measurements, and 
measurements of the amount of deposition near the pad, were made after the first launch 
studied, STS 41 D. For the second launch, STS 51 A, only HClgas concentrations were 
measured. Deposition measurements for the 41 D launch will be discussed first, followed 
by a description of the gas and aerosol measurements. 

The same oil-cup technique developed for the model tests was used to sample the 
acid deposition on and around the launch pad at KSC. For safety reasons personnel must 
clear the pad area before fueling of the external tank begins on the night prior to the 
launch, and they cannot return until the pad is cleared several hours after lift-off. Thus 
the oil cups were deployed and left uncapped from prior to 8 p.m. the night before and 
picked up beginning shortly after noon on the day of the launch. Figures 12 and 13 show 
the deployment array. On Figure 12 the sites are annotated to indicate the appearance 
of the oil and the depth of the solid material collected; on Figure 13 they are annotated 
with acid normality (when the sample was large enough to obtain a measurement) and 
liquid depth. Both clear and greenish deposition drops were noted in the samples. The 
milky appearance of the oil in some cups was probably due to collection of fine liquid 
spray from the launch; the milkiness was greatest at locations where one might expect 
the most intense spray. In some samples the oil remained clear, so the possibility that 
the milkiness was caused by dew fall can be discounted in this case. 

The HC1 gas measurement effort was undertaken as a joint project sponsored by 
NASA and the Air Force Engineering and Services Center, Environmental Sciences 
Branch, Tyndall Air Force Base, Florida. The work was performed by the Arnold En- 
gineering Development Center (AEDC), Air Force Systems Command. Michael G. Scott 
and Charles W. Pender, Jr., of Sverdrup Technology, Inc. (operating contractor for the 
propulsion test facilities at AEDC) were the principal investigators. Lt. Floyd Wiseman, 
AFESC/RDVS, was the Project Officer for the Engineering Services Center and Capt. 
Frank Tanji, AEDC/DOTR, was Project Manager for AEDC. 

A long path infrared absorption technique using a Fourier transform spectrometer 
(FTS) was selected for this study because the launch site is on the sea coast and the 
background chloride concentrations may be high at times. Thus a gas measurement 
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Figure 11. Deposition acidity pattern for the August 26, 1983, VLS configuration model test. For 
this test, the model was anchored in an elevated position to simulate lift-off conditions. 



Figure 12. Deposition pattern found after the STS 41 D launch. 

The line of site for the Fourier transform spectrom- 
eter gaseous 1IC1 measurement is also indicated (FTS 
to blackbody). 
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Figure 13. Liquid deposition depth and acidity pattern for the 
STS 41 D launch. Depth is indicated in centimeters. 



technique which is highly specific for HC1 was required. In this technique, an infrared 
light source (a high temperature blackbody) is set up at one location and viewed with a 
teiescope-FTS system from several hundred meters away. The spectrometer measures, 
as a function of wavelength, the absorption of the infrared energy caused by the gases 
along the path from the blackbody. Since HC1 and other gases absorb at specific, known, 
wavelengths in the infrared, the presence of HC1 can be definitely determined and the 
quantity measured from the ratio of the strength of the absorption lines to the back- 
ground envelope. Laboratory calibration showed that the detection threshold for the 
system was about 0.5 ppm by volume. 

At KSC the exhaust from both SRBs is ejected from a single horizontal trench 
which opens to the north onto a flat grass-covered area. The grass extends approximate- 
ly 400 m from the launch pad to a lagoon and brush-covered area beyond the launch 
complex perimeter. The bulk of the wet acidic deposition from each launch, falls onto 
this grass-covered area and the lagoon and brush beyond, making this the primary source 
of HC1 revolatilization. Measurements for both the 41D and 51A launches were made 
directly over the grass-covered area within the perimeter fence. The absorption was 
measured along a 500-m path oriented east-west, approximately 1.2 m above the surface. 
The blackbody was east of the HC1 source region; the FTS was on the west side. From 55 
to 60 percent of the path was directly over the HC1 source (Fig. 12). 

The FTS used in this study was a Block Model RS197 configured for maximum 
sensitivity in the spectral region of interest, 3000 to 2700 cm . The system used a 
germanium beam splitter on a potassium bromide substrate and an indium antimonide 
detector cooled to 77 K. Special precautions in the mounting and housing of the instru- 
ment were observed to enable it to withstand the severe acoustic vibration and corrosive 
HC1 gas environments produced by the shuttle launch. At regular intervals, the system 
would collect 124 interferograms in a 2-min period, average the digitized data, and store 
the results on magnetic tape. The averaged interferograms were transformed into the 
spectral domain and the amount of HC1 absorption determined at a later time. Sample 
absorption spectra from the laboratory calibration and the field data are illustrated in 
Appendix I. The interval between data collection was 20 min for the 41 D launch; soft- 
ware improvements allowed this to be reduced to 10 min for 51A. For additional details 
on the system design, calibration, and operation, the reader is referred to Appendix I. 

Hydrogen chloride revolatilization was measured after the 41D and 51 A launches 
as shown in Figure 14. In Table 4, meteorological data and other relevant information 
concerning these launches are summarized. The HC1 gas concentrations plotted as 
functions of time in this figure are computed assuming the gas is evenly spread through- 
out the 500-m path between the blackbody radiation source and the FTS receiver. In 
actuality, the concentration is expected to be somewhat greater over the source and less 
elsewhere. The precision of the 41 D data is ±19 percent, ±15 percent for 51 A. The 
difference is due to changes in the blackbody radiation source made between launches 
which resulted in a factor of 4 improvement in the signal-to-noise ratio. 

For both launches, the 1IC1 concentration reached a peak just over an hour after 
launch, remained high for an hour or so, then slowly decayed. The measurements were 
continued for 3 days following each launch. After 41 D, no HC1 was detectable after the 
7 -hr period illustrated by Figure 14. This was not the case after 51 A. Trace amounts, 
less than a part per million, were detected between 8 a.m. and 2 p.m. local time on the 
day following the launch, and again between 10 a.m. and 1 p.m. on the second day after 
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SHUTTLE LAUNCH 41 D & 51 A 





the launch. The improvement in signal-to-noise ratio may account for the difference, 
but the initial source strength was also stronger after 51A. 


TABLE 4. ENVIRONMENTAL CONDITIONS FOR 41D AND 51A LAUNCHES 


Condition 

STS 41 D 

STS 51 A 

Launch Date 

Aug. 30, 1984 

Nov. 8, 1984 

Launch Time (local) 

08:42 EDT 

07:15 EST 

Temperature (Celsius) 

26 

20 

Relative Humidity (percent) 

81 

59 

Wind Speed at 60 ft. (m s - *) 

5 

7 

Wind Direction (degrees) 

106 

24 


NOTE: The ground was dry during the 41 D launch except for possible due from 
the night before. After the 51A launch, the ground was observed to be 
quite wet with standing water in some locations, apparently from 
thunderstorms in the area during the night. 


Since the measurement method discussed above only detects gaseous HC1, the 
possibility that HC1 vaporizes from the ground and then condenses to an aerosol phase 
must also be considered. In the aerosol phase, HC1 in concentrations of a few parts per 
million should be readily detected by standard aerosol counting and sizing techniques, 
assuming there is no undue interference from local anthropogenic sources. For example, 

1 ppm HC1 by volume is equivalent to about 1.5 mg HC1 per cubic meter. On the 
coast at ground level, sea salt aerosol concentrations of 1 to 15 yg m”'* for wind speeds 
up to 8 m s -1 may be expected (13). If the relative humidity is high, the background 
aerosol mass will increase by deliquescence, but it must exceed 98 percent relative 
humidity before the background would match the HC1 mass concentration. Thus the HC1 
aerosol should, as a minimum, double the background levels if the concentration is 1 ppm 
or more in aerosol form. 

Aerosol number concentration and size distributions were measured after the 41 D 
launch using a Gardner counter to provide the total Aitken nucleus count (particles 
between 0.01 and 0.2 y m) and a Climet optical particle counter to give the number and 
size distribution for particles between 0.3 and 10 y m. The raw pulse output from the 
Climet was fed to a multichannel analyzer so that the entire size spectrum could be 
determined for each sample. The spectrum accumulated by the multichannel analyzer 
was dumped to paper tape at regular intervals. The system was housed in a van so that it 
could be moved alternately upwind and downwind of the source area. 

It was intended that the Climet system would operate automatically, beginning 
the night before, right through the launch and post-launch period. Unfortunately, an 
equipment malfunction prevented this from happening so data were obtained only while 
the instrument was manually attended, beginning 1:49 p.m., after access to the pad was 
allowed. (The launch occurred at 08:42 EDT.) Eight 70.8-liter air samples, each drawn 
over a 10-min period, were taken downwind from the source area (at the FTS instrument 
location) between 1:49 and 3:53 p.m. Then between 4:30 and 5:30 p.m., five samples 
were taken upwind, at the camera site on the east side of the pad. Finally, four addi- 
tional samples were collected downwind at the FTS site between 5:40 and 6:40 p.m. With 
the exception of one sample drawn just after 2 p.m., there was no systematic difference 
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between the up- and downwind samples; the difference between the average counts was 
less than the standard deviation. The count in the sample taken just after 2 p.m. was 2.8 
times the mean, and the following sample also showed a high count, but traffic in the 
area was high at the time so the high counts are probably due to road dust. 

The question of HCl/f^O aerosol formation has been treated analytically by Rhein 
(14) and, more extensively, by P.V.N. Nair et al. (15). Their work indicates that at 20 °C 
the threshold for HC1 aerosol condensation is 10 ppm at 80 percent relative humidity and 
1 ppm for relative humidities over 91 percent. Thus examining the gas concentration 
data for the times of the aerosol measurements, aerosol formation would not be expected 
under these conditions, in agreement with the aerosol measurements. 

IV. HC1 REVOLATILIZATION ANALYSIS 

The qualitative features of the HC1 revolatilization process as illustrated by the 
51 A and 41 D data (Fig. 14) may be explained quite simply. As discussed above, the 
launch leaves the pad area and the immediate vicinity on the north (SRB exhaust) side of 
the pad covered with drops or small pools of hydrochloric acid. Initially the acid is fairly 
dilute, about 2 N or less. For these low concentrations, the equilibrium vapor pressure of 
HC1 is more than five orders of magnitude less than that of water over the solution. 

Thus the initial acid vapor concentration in the air is relatively low. However, it in- 
creases rapidly as the drops, small ones first, lose water by evaporation and increase in 
acidity. When evaporation has proceeded to the point that less than a quarter of the 
mass of a drop remains, the concentration reaches a value on the order of 11 N and the 
acid and water vapor pressures are equal. Thus most of the acid is released in the final 
stage of evaporation of each drop. Therefore, the ambient HC1 vapor concentration rises 
to a peak as the majority of drops evaporate, small ones first. Then it falls slowly, 
fueled only by the much slower evaporation of the largest drops, small pools, and acid in 
the surface soil layer. 

To develop a quantitative treatment of this process we note first the great sim- 
ilarity to the problem of evaporation of water or other pure liquid from a pool or field. 
These problems have been treated extensively in the literature because of their great 
importance to hydrology, agriculture, and the study of chemical spills. A full treatment 
is quite complex, since the evaporation rate depends on many variables including wind 
speed, ambient humidity, solar insolation and cloud cover, atmospheric stability and 
turbulence, terrain features and surface roughness, vegetation cover, etc. The HC1 
revolatilization is more complex, however, because it involves the interdependent evapo- 
ration of two substances, acid and water. 

In this study, an attempt is made only to elucidate the basic nature of the re- 
volatilization process and the measurements already discussed. Development of a 
detailed revolatilization model is beyond the current scope of the project. Instead, a 
highly simplified treatment of the key aspects of the problem has been developed which 
will serve to explain the basic physics. This treatment is expressed as a simple numerical 
model listed in Appendix II. The essential aspects of the model are as follows: 

1. The source area is assumed to be covered initially with a Gaussian distribution 
of hemispherical drops, all of the same acid concentration. Initial parameters to be 
specified: mean radius, total liquid volume per square meter, standard deviation of size 
distribution, and acid content (weight percent). 
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2. The rate at which pure water would evaporate from the surface is computed 
using a simple empirical formula defined in the classic work by Penman (16). The evapo- 
ration rate is expressed as a function of the following variables: wind speed, ambient 
humidity, air temperature, and surface temperature. The surface temperature is typical- 
ly not measured and it is fairly difficult to calculate since it is a function of a number of 
parameters. Specifying it as an input is one of the major simplifications of this treat- 
ment. 


3. Beginning with the initial conditions, the evaporation of acid and water is 
calculated as a function of time. The evaporation from each drop size class is de- 
termined by the number of drops and the microphysical parameters, drop radius, acid 
concentration, surface temperature, and an estimate of the near-field acid and water 
vapor concentrations from the prior time step. At each time step, the total amount of 
material evaporated is normalized via the ratio of latent heats to the amount of evapora- 
tion (water) given by the Penman equation. 

4. As the evaporation process proceeds and the surface begins to dry, the total 
amount of evaporation given by the Penman equation is adjusted down in proportion to 
the ratio of wet to dry surface. 

The model outputs the HC1 source strength per square meter of surface as a 
function of time. The mixing process into the atmosphere is not modeled. For this 
reason, and because ground temperatures were not measured following the launches for 
which data were obtained, a direct comparison between model output and measurements 
is not possible. However, comparing Figure 15, which displays the model output for cases 
approximating the 41 D and 51 A post-launch conditions, to Figure 14, it is clear that at 
least the basic form of the functional dependence with time is correct. If the HC1 from 
a strip 100 m long by 1 m wide (aligned with the wind) is assumed to mix into a volume of 
air 2 m high by 1 m wide by U x Dtime long (U = wind speed and Dtime = the time 
interval), it is found that the model output is the right order of magnitude to explain the 
observed results. Examination of Table 5 and Figures 16 through 22, which illustrate the 
dependence of the HC1 source on the various parameters, shows that the model results 
are quite sensitive to these variables. Thus considerable additional work is required to 
complete validation of this analysis. It is presented here to indicate the types of depen- 
dence the HC1 revolatilization is expected to show on the various parameters and to 
serve as a basis for future modeling efforts. 

It is anticipated that a complete model of the revolatilization process could be 
developed quite readily based results is desired, it could be obtained by replacing the 
empirical Penman equation with a simple model for evaporation from a pool such as the 
modified Ille and Springer model as discussed by Kunkel(17); coupling the results (the 
source term) to a Gaussian diffusion model like that discussed in Refs. (18) through (20). 
Additional work to improve the numerical methods and reduce the computation time 
required by the current model would be well worthwhile. 

V. SUMMARY 

Whenever large solid rocket motors which produce hydrogen chloride as an exhaust 
product are launched or test fired from a facility which causes the exhaust to mix with 
large volume water sprays (more than can be evaporated by the exhaust heat), an ap- 
preciable acid deposition in the near field is going to result. Some details of the produc- 
tion mechanism and, more importantly, the scaling are not yet fully worked out, so 
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Figure 15. Analytical model IIC1 source strength: (1) STS 41 D conditions; (2) STS 51 A conditions. 




TABLE 5. CROSS TABLE OF PARAMETER VARIATIONS FOR 


Case 

Designation 


BAS47L 

(Baseline) 

U-VAR 1 
U-VAR 2 
U-VAR 3 
D-VAR 1 
D-VAR 2 
D-VAR 3 
T-VAR 1 
T-VAR 2 
T-VAR 3 
A-VAR 1 
A-VAR 2 
A-VAR 3 
R-VAR 1 
R-VAR 2 
R-VAR 3 
V-VAR 1 
V-VAR 2 
V-VAR 3 
SGVAR 1 
SGVAR 2 
SGVAR 3 



**No wind. 


NOTE: The baseline case is the best fit to the STS 41 D post-launch 

data. The actual time variations of air temperature, dew 
point, and wind speed were approximated by linear curve fits. 
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TIME (minutes) 

Analytical model HC1 source strength results illustrating dependence on area factor for 
evaporation reduction: (1) BAS47L; (2) A-VAR 1, 1; (3) A-VAR 2, 3; (4) A-VAR 3, 4. 
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Figure 20. Analytical model 1 1 Cl source strength results illustrating dependence on drop radius: (1) 
BAS47L; (2) R-VAR 1, 0.05 cm; (3) R-VAR 2, 0.15 cm; (4) R-VAR 3, 0.20 cm. 
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TIME (minutes) 

Figure 22. Analytical model IICl source strength results illustrating dependence on the width of 
the drop size distribution, sigma: (1) SGVAR 1, 0.02 cm; (2) SGVAR 2, 0.04 cm; (3) 
SGVAR 3, 0.06 cm; (4) BAS47L. 



quantitative predictions of the amount and location of the deposition can only be esti- 
mated for new systems. However, the important parameters have been identified and 
the estimates should be adequate for most engineering applications, especially since the 
wind and other weather effects may be expected to generate considerable variability in 
the initial deposition pattern and subsequent effects. 

The life cycle of the acid deposition, once it occurs, is well established by the 
field and model scale measurements and the analysis presented here. The initial deposit 
is moderately acidic, of order 2 N for space shuttle launches at KSC. As the material 
evaporates, the water vaporizes much more rapidly than the HC1 so the acid 
concentrates until the water and HC1 vapor pressures are equal, at about 11 N. If the 
evaporation potential is high (warm temperatures, low humidity, and moderate to high 
wind speeds), the evaporation will proceed rapidly enough that corrosion damage from 
direct contact with the liquid will not be immediately evident except on the most 
susceptible exterior surfaces. The dominant effect in this case will most often be the 
aluminum oxide particulate deposition which is scavenged from the exhaust along with 
the HC1. The acid has the effect of increasing the bonding between the aluminum oxide 
and structure surfaces, so the surface ends up coated by a powder which is difficult to 
remove except by direct mechanical scrubbing. At launch or test sites where this 
coating is expected to be a problem, the addition of chemical additives to the facility 
water to reduce this bonding should be investigated. 

Under meteorological conditions when the evaporation potential is low, experience 
has shown that the concentrated deposition remains on painted and metal surfaces long 
enough to cause immediate corrosion damage, spotting on automobile chrome for 
example, and burn spots on vegetation. The residual aluminum oxide powder is still a 
problem under these circumstances, although the drying may be slow enough that the 
timely application of sprinkler systems and washdown hoses may alleviate the situation 
on the portions of the facilities that can be reached before evaporation is complete. 

For Vandenberg SLC-6 and similar facilities where extensive computer and elec- 
tronic equipment is located in close proximity to the launch pad, the most serious prob- 
lem associated with acid residue from a firing is not likely to be the liquid deposition 
itself, but the HC1 gas which evolves as the liquid evaporates. Equipment of this type is 
often very sensitive to corrosion damage from gas concentrations in the 10 to 100 parts 
per billion level. Exposures of 8 to 10 hours may render a computer system inoperable. 
The measurements reported here confirm that HC1 concentrations above this level may 
be expected intermittently at the launch site for at least 2 days following a firing. They 
also verify that the concentration can exceed 5 ppm for brief periods in the first few 
hours; 5 ppm is the threshold limit value for workers; 1 ppm is the public exposure limit. 
Dangerous levels in low, enclosed, or partially enclosed structures are also a possibility. 
Thus the safety aspects of the HC1 evaporation must not be ignored, although fairly 
straightforward precautions should be adequate for most situations. 
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PREFACE 

The work reported herein was conducted by the Arnold Engineering Development Center 
(AEDC), Air Force Systems Command (AFSC), at the request of Air Force Engineering 
Services Center (AFESC/RDVS), Tyndall Air Force Base, Florida, and the Atmospheric 
Effects Branch, Marshall Space Flight Center, NASA H-73999B. The reported measurements 
were conducted and results obtained by Sverdrup Technology, Inc., AEDC Group, operating 
contractor for the propulsion test facilities at the AEDC, AFSC, Arnold Air Force Station, 
Tennessee, under Project Number DB07EW. The Air Force Project Managers were Captain 
Frank Tanji and Captain Bradley Biehn, AEDC/DOTR. The data analysis was completed 
on May 15, 1985, and the manuscript was submitted for publication on July 1, 1985. 
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1.0 INTRODUCTION 

Large amounts of gaseous hydrogen chloride (HC1) are one of the combustion by-products 
generated during the launch of the Space Transportation System (Shuttle) by the Shuttle’s 
solid-propellant rocket boosters (SRB’s). Significant quantities of HC1 are present in the plume 
cloud that lingers near the launch structure. The HC1 gas remains entrained in the atmosphere 
near the launch facility, is dispersed by wind, and is deposited on surfaces, such as grass 
or ground. (See Ref. 1 for deposition mechanisms.) The HC1 gas that remains in the vicinity 
of the launch structure, from either initial release or secondary ground release, is important 
because it corrodes metals and electronic equipment and because it is a health hazard. The 
problem is currently restricted to the Kennedy Space Center (KSC) launch area but will also 
be present at Vandenberg when the Shuttle Western Test Range launch facility is activated. 
The problem may be exacerbated at Vandenberg because of differences in the launch mount, 
sound-suppression water systems, and the proximity of the launch facilities. 

The objective of this project was to develop and apply a nonobtrusive absorption technique 
employing a Fourier Transform Spectrometer (FTS) to monitor the concentration of gaseous 
HC1 present near the KSC launch pad for several days after a Shuttle launch. Because of 
safety considerations and the requirement to collect data immediately after a launch, the 
measurement technique was automated to permit untended operation. The technique was 
applied for launches 41D and 51A, Pad 39A (Fig. 1), which occurred on August 30, 1984 
and November 8, 1984, respectively, at KSC. 

To accomplish the objective, a technique based on the absorption of infrared (IR) radiation 
by the diatomic HC1 molecules was developed. The characteristic HC1 lines are attributable 
to absorption of radiation associated with vibrational and rotational transitions within the 
HC1 molecule. A simple model (Ref. 2) to describe this absorption process is a molecule in 
which the individual atoms, held together by chemical bonds, are in vibratory motion along 
these bonds, while the entire molecule is rotating. The HC1 is in a state of vibratory motion 
brought about by the alternate stretching and contracting of the chemical bond as the hydrogen 
and chlorine atoms move away and toward each other, respectively. This vibratory motion 
is superimposed on a rotation of the molecule about an axis perpendicular to the chemical 
bond. When IR radiation of the proper frequency (i.e. energy) impinges on the molecule 
and is absorbed, the vibration and/or rotation states are changed. These changes must satisfy 
certain selection rules that give rise to discrete absorption lines. These lines are labeled according 
to their frequency. The lines with frequencies greater than that of the band center are said 
to be in the R branch, whereas those with frequencies less than that of the band center are 
in the P branch (Fig. 2). The individual lines making up the P and R branches are identified 
as P(l), P(2), etc. as the frequency of the lines moves away from the band center. 
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2.0 TECHNIQUE 

A nonobtrusive IR absorption technique was used for monitoring gaseous HCI 
concentrations over long paths near the Shuttle launch pad (Fig. 1). The technique involved 
use of an IR spectrometer located approximately 500 m from an IR source. An IR spectrum 
was acquired at regular intervals (20 min for 41D, 10 min for 51A) for severed days following 
the Shuttle launches. The spectra were then inspected to find the magnitude of absorption 
because of the presence of gaseous HCI. The concentration of HCI was then calculated using 
the absorption coefficient measured in the laboratory. The spectral line used to determine 
the amount of absorption is in the band centered at 2,886 cm -1 . This technique is not 
dependent on absolute magnitudes of the spectra; therefore, the intensity of the source is 
not relevant to the measurement, and only relative intensities will be presented in this report. 

Naturally occurring HCI contains two isotopes of chlorine. They are 35 C1 and 37 C1, in 
the ratio of 3 to 1 , respectively. The IR absorption lines of H 35 C1 and H 37 C1 are separate, 
but very close. A spectrum of these isotopic species is shown in Fig. 2b. 

The absorption attributable to HCI gas is well understood and documented. The HCI 
absorption lines are described by (Ref. 3) 

I(\) = I 0 (X) exp { -SP 7 o L/[(X - X 0 )2 + ^ 


where 

I Q (X) is the intensity of the incident radiation 

I (X) is the intensity of the transmitted radiation 

S is the line strength of the absorption line (cm -2 atm -1 ) 

P is the partial pressure of the absorbing species (atm) 

7o is the pressure-broadened half-width at half-maximum of the HCI line (cm -1 ) 

Xq is the line center of the absorption line (cm -1 ) 

L is path length 

Figure 3 shows the P(l) lines calculated using parameters contained in Ref. 4. The figure 
shows the theoretical absorption assuming 6 ppm HCI over a 500-m path. Because of limited 
resolution and the apodization function, the Fourier transform spectrometer (FTS) that 
produced the actual spectra influences the spectra shape. The convolution of the instrument 
line shape (ILS), using 0.5-cm - 1 resolution and a triangular apodization function, produces 
the P(l) line shown in Fig. 3b. As can be seen by comparing Figs. 3a and b, the ILS shortens 
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and broadens the spectral lines. (Note: The scales of these plots are not the same.) The effect 
is detrimental in that the resultant spectrum is not an accurate portrayal of the true line shape. 
However, it is still possible to use the resultant spectra since the effect can be calibrated. 
The calibration will be discussed later in this report. 

A complicating factor was the presence of additional absorption lines in the spectra obtained 
in the field. The effect of neighboring absorption lines is to give the impression of a greater 
amount of gaseous HC1. The spectral region of interest included absorption lines attributable 
to other molecules present in the atmosphere such as H 2 0 and. Hydrogen Deuterium Oxide 
(HDO) (Ref. 5). Therefore, the spectra were obtained at 0.5-cm -1 resolution to separate 
the HC1 lines from the lines of interference. In Fig. 4a the HC1 doublets can be easily discerned 
when compared to a background spectrum with no HC1 present. Fig. 4b. Figure 4a is a portion 
of a typical spectrum acquired during the 51 A launch. This figure shows only the lines within 
the P branch. The R branch is unusable because of domination by water vapor bands at 
2.7 and 3.2 jim. The P(5) line was used to determine the HC1 concentrations because it had 
the least interference. 


3.0 APPARATUS 

The Shuttle launches took place at KSC Pad 39A. Figure 1 shows the major features 
of the launch pad and the AEDC equipment location. The test apparatus consisted of an 
FTS and a collimated light source with associated control and data acquisition instrumentation. 
The sample path length was 500 m. The distance to the launch structure was approximately 
350 m. The equipment was located so that the plume deflected by the flame trench passed 
through the center of the sample path. 

The FTS used was a Block Model RS197 field-rated instrument (Refs. 6 and 7). It was 
configured to have the maximum sensitivity possible in the spectral region of interest, 3,000 
to 2,700 cm -1 . This was accomplished by using a germanium beamsplitter on a potassium 
bromide substrate and an indium antimonide detector cooled to 77 K. The detector dewar 
had a maximum hold time of 6 hr, so an automatic fill system was developed to replenish 
the liquid nitrogen (LN 2 ) from a large auxiliary reservoir. The FTS was rigidly mounted to 
maintain optical alignment with the source during the large vibrations associated with the 
Shuttle launch. The instrument was enclosed in a Plexiglas® container that was kept at a 
slight positive pressure with a gaseous nitrogen (N 2 ) purge to reduce the possibility of HC1 
corrosion (Fig. 5b). The IR radiation entered the FTS through a hole in the Plexiglas container 
that matched the entrance aperture of the FTS. An O-ring seal was made between the FTS 
and the container. In this way, the FTS was in the same configuration in the field test as 
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in the calibration. The only difference in the optical paths was that the quartz sample cell 
was fitted with IR transmisive quartz windows on each end. 

The source units used in launches 41 D and 51 A were different. Modifications were made 
to the equipment for launch 5 1 A because of insight gained during the 41D launch. The original 
configuration (See Fig. 6a) used in the 41D launch consisted of a Barnes Model 112017 
blackbody with collimator and an Optronics Model 100U, 1,000-w tungsten halogen lamp 
with a 25.4-cm Cassegrain telescope used as a collimator. The use of multiple sources was 
required to provide adequate intensity since the FTS was located 500 m distant and also as 
insurance against optical misalignment from the shock of the booster ignition. In the second 
configuration (See Fig. 6b), for the 51A launch, three, 30.48-cm spherical mirrors with 1,000-w 
tungsten halogen lamps at their focal points were used to provide an intensity increase of 
400 percent. In each case the sources were securely mounted and covered with a hood that 
provided protection from inclement weather and SRB plume debris. The sources were mounted 
so that the optical path between the FTS and sources was approximately 1 .2 m above the 
ground (Fig. 6b). 

The digital data acquisition system collected a set of 124 interferograms (the raw data) 
during each sampling period of approximately 2 min. Each set was averaged and stored on 
magnetic tape as one data point. A spectrum is obtained by taking the Fourier transform 
of the interferogram (Ref. 6). For launch 41D, data points were collected every 20 min. This 
time interval was reduced to 10 min for 51 A after improvements were made to the software 
that reduced the processing time from 17 to 8 min. The data were not converted to the spectral 
domain until later in order to permit the collection of additional data during the postlaunch 
period. 


4.0 CALIBRATION AND ANALYSIS 

The absorption technique involved the comparison of spectral data acquired in the field 
with similar data acquired in the laboratory. As explained in Section 2.0, calibration data 
generated in the laboratory were used in lieu of published HC1 absorption coefficients to 
negate the influence of the spectrometer instrument line shape. 

In the laboratory (Fig. 5a), a previously evacuated quartz sample cell (2.54 diam by 25.4 
cm) was filled with mixtures of HC1 and dry N 2 with known concentrations simulating the 
absorption expected following the launch. The concentration of HC1 in the sample cell is 
stated in terms of partial pressures (x) with units of psia, and the concentration (y) in the 
field is expressed as parts per million (ppm). The relationship is x/14.7 = y/1 ,000,000. Specific 
sample cell concentrations of HC1 were achieved using a partial pressure technique in which 
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the cell was filled in several stages. Each stage consisted of partially evacuating the cell followed 
by bringing it back up to 14.7 psia with dry N 2 . This process was continued until the desired 
partial pressure of HC1 was reached. To simulate the field conditions, the concentrations 
used in the laboratory were greater by an amount dictated by Beer’s law. A simple form 
of Beer’s law. 


ln(I/I 0 ) = -ABC 


where 


I = transmitted radiation 
I 0 = incident radiation 

A = absorptivity (constant at a given wavelength) 

B = path length 

C = concentration of the absorbing medium 

states that the ratio of I/I 0 will remain constant if C is varied inversely to B. For example, 
I/I c is equivalent for a gas at either 1 psia over a 1-m path or for 10 psia over 0.1 m. 
Therefore, in the laboratory the concentrations needed to be greater by a factor of 1,968 
(the ratio of the range in the field, 500 m, to the cell length, 25.4 cm). 

A set of nine spectra corresponding to nine calibration pressures were obtained. The x 
axis of these curves is in terms of path-averaged concentrations over 500 m. Calibration curves 
were obtained by plotting ln(I 0 /I) versus x for six P lines. The plots are shown in Figs. 7a 
through f. The partial pressures of HC1, in psia, used in the test cell were 0.018, 0.036, 0.053, 
0.080, 0.120, 0.181, 0.272, 0.544, and 1.633. The conversion relationship to the equivalent 
field concentration is given by y (ppm) = 34.6 x(psia). 

The analysis of the KSC spectra consisted of measuring the difference between the baseline 
and the peak depth of the six P absorption lines. An important feature of this technique 
is that atmospheric haze does not invalidate the measurement. Atmospheric scattering 
attributable to haze does not affect the procedure since both I and I c (Fig. 2b) are influenced 
by the same amount, leaving the ratio I/I 0 constant. This feature negates the necessity of 
performing an instrument or source calibration. 

The error associated with the resultant measurements of the 41D launch is estimated to 
be ± 19 percent (root sum of the squares). This error is caused by an error of ± 10 percent 
because of pressure transducer uncertainty, ± 10-percent gas-handling technique, ± 10 percent 
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in FTS measurement uncertainty, and ±8 percent in data handling (digitization and 
processing). The error associated with the resultant measurements of the 51 A launch is 
estimated to be ±15 percent (root sum of the squares). This reduction in error is caused 
by the use of the alternate source unit with a higher intensity during the second measurement. 
The FTS measurement uncertainty was reduced to ±3 percent, and the data handling 
uncertainty was reduced to ±2 percent. 

5.0 RESULTS 

The application of the IR absorption technique following the 41D launch showed the 
presence of average concentrations of gaseous HC1 as great as 4 ppm ( ± 19 percent). Figure 
4a is a representative spectrum taken postlaunch. The HC1 lines are clearly apparent among 
other atmospheric absorption lines. Figure 8 gives the concentration of gaseous HC1 as yielded 
by examination of five H 35 C1 absorption lines. The spread in the results is attributed to the 
presence of absorption lines of other molecules. Note that in Fig. 8 the concentration 
determined using the P(5) line is less than concentrations determined using the other lines. 
In the case of most of the HC1 absorption lines, their frequencies are nearly the same as 
other atmospheric constituents, such as HDO and H 2 O (Ref. 5). Therefore, the resultant 
absorption is greater than that caused only by the H Cl. Since it is difficult to determine how 
much of the absorption is caused by the other molecules, the most isolated line, the P(5) 
line, was used to determine the HC1 concentration. 

The concentration peaked approximately 100 min postlaunch, remained high for nearly 
1 hr, and then began to diminish. There unfortunately exists a gap in the data because of 
a period in which instrumentation was realigned following the loss of one of the sources 
during the high vibration experienced during the launch. The instrumentation was operated 
for three days following the launch, and HC1 was detected only in the 6-hr interval following 
the launch. 

Data obtained from the 51 A launch are better in two ways. Data were acquired more 
frequently and with an improvement in the signal-to-noise ratio by nearly a factor of four. 
During the 51 A launch, higher HC1 concentrations were measured than were detected during 
the 41D launch. Figure 9 is a summary of six concentration curves from the P-branch lines. 
The curves have a significant scatter caused by absorption attributed to other molecules present 
in the sample path. The most isolated line again is the P(5) line. Figure 10 shows the best 
estimate of the HC1 concentration. The HC1 concentration peaked at about 9 ppm ± 15 percent 
approximately 1 hr postlaunch, fluctuated for about an hour, and then decayed. In the two 
days following the launch, HC1 was detected in minute quantities after sunrise until early 
afternoon. Long-term, low levels of HC1 were visible in the 51 A launch either because of 
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greater revolatilization attributable to atmospheric and ground conditions or better signal- 
to-noise ratio. Atmospheric conditions during the measurement periods of both launches 
are available (Refs. 8 and 9). 


6.0 SUMMARY 

The objectives of the project were completely fulfilled. An unobtrusive absorption 
technique to monitor concentrations of gaseous hydrogen chloride was developed and then 
applied on two space shuttle launches at the Kennedy Space Center. The technique developed 
exploited the intrinsic property of gaseous HC1 to absorb IR radiation in specific, narrow, 
spectral lines. This IR absorption is measured by a Fourier Transform Spectrometer that 
features accuracy and precision without a calibrated source. The minimum detectable 
concentration of the instrument is 0.10 ppm with an error of 15 percent (obtained by taking 
the square root of the sum of the squares of all possible errors or uncertainties in the 
measurement system including calibration). The application portion of this project enabled 
refinement of field use of the measurement technique. Increased source radiance and larger 
source collimating optics improved the signal-to-noise ratio by a factor of four. Improvement 
in data processing software decreased uncertainty by 4 percent. 


Measured quantities of gaseous HC1 following launches of both 41 D and 51 A displayed 
similar trends. The concentration after each launch first increased, then decreased, then 
increased to a maximum approximately 100 sec after the launch. The concentration decreased 
gradually to less than 1 ppm approximately 10 hr after each launch. The peak concentration 
measured was 4 and 9 ppm for 41D and 51A, respectively. For two days after the 51 A launch 
only, detectable levels were recorded for a period of approximately 6 hr following sunrise. 
Atmospheric and ground conditions were different for each launch (Refs. 8 and 9). Additional 
considerations in interpreting the concentration levels are that the data are both path and 
time averaged, and that the ground scar from the plume makes up approximately 60 percent 
of the 500-m sample path length. These considerations imply that the concentration along 
the path is not uniform. Additionally only gaseous HC1 will be detected by this technique. 

7.0 RECOMMENDATIONS 

Improvements to the instrument developed on this project can be made in two areas. 
The signal can be increased by enlarging the collection optics on the FTS, and the noise can 
be reduced by increasing the data averaging period. Some improvement in calibration accuracy 
could be realized by using gas-mixture ratios prepared by the National Bureau of Standards. 

The technique developed by this project could use computed tomography with multiple 
paths to generate a two-dimensional map of HC1 concentration levels to determine areas of 
highest concentration and change over time. This technique would enable determination of 
safe areas around solid-propellant rocket motor launch complexes where people, wildlife, 
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or machinery may be adversely impacted. Additional analysis (not within the scope of this 
project) should be conducted to determine the impact of atmospheric and ground conditions, 
aerosol entrapment of HC1, and flame trench cooling water runoff and holding pond 
contributions to gaseous HCI concentration levels. 
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HC1 Absorption Sepctrum 



a. Covering 3,200 to 2,200 cm-i 
Figure 2. Sample spectrum showing HC1 fundamental 
absorption band. 
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HC1 Absorption Spectrum 



b. Showing only the P branch and how I c and I are determined 
Figure 2. Concluded. 
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a. Typical spectrum obtained after the 51A launch. 
Figure 4. 51A launch spectra. 
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b. Background spectrum obtained before the 51A launch 
Figure 4. Concluded. 
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Time Following Launch in Minutes 


Figure 8. Summary of 41D HC1 concentrations measured using several 
P-branch absorption lines. 
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Figure 9. Summary of 51 A HC1 concentrations measured using several 
P-branch absorption lines. 
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Figure 10. 51A HC1 concentration using P(5) line for 10 hr 
following launch. 
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Analysis Program for Estimation of HC1 
Revolatilization Source Strength 

REM HCL REVOLITIL3Z ATION PROGRAM: VERSION "HCL41L" 

REM INCORPORATES REDUCTION OF EVAPO IN PROPORTION TO SUMAREA# 

REM INPUT MODIFIED TO MATCH 41 D LAUNCH, POST LAUNCH CONDITIONS • 
REM USING LINEAR CURVE FIT FOR Ts, Tdew, AND U* 

REM WRITTEN BY JEFFREY ANDERSON, OCTOER 16, 1986 
OPTION BASE 1 
INTEGER Loop 

DIM N(40) ,A(40) ,C(40) ,V(40) ,Mga(40) ,M(40) ,Mtg(600) ,Mw( 40 ) , Pow ( 40 ) 
DIM Case$ [60] 

INPUT "IS INITIALIZED DATA TAPE IN DRIVE :T14? (Y/N)”,A$ 

IF A$=" Y" THEN 110 
GOTO 80 

REM N ( 1 ) = number of drops of radius A(I) 

REM A ( I ) = DROP RADIUS IN CM, DROP IS ASSUMED HEMISPHERICAL 
REM C(I) = HCL CONCENTRATION IN WEIGHT PERCENT 
REM V(I) = VOLUME OF DROP IN CUBIC CM 

REM Mga(I) = MASS OF HCL GAS EVAPORATED FROM Ith DROP (GRAMS) 
REM M ( I ) = MOLARITY OF DROP 

REM Mtg(T) = TOTAL HCL (grams) EVAP IN THE TIME STEP 

REM Mw( I ) = MASS OF WATER EVAPORATED FROM ALL RADIUS I DROPS 

REM 

REM 

REM DEFINE CONSTANTS AND INITIAL CONDITIONS 

REM 

Dif f w=* 25 ! DIFFUSION COEF OF WATER VAPOR, CM^2 /S 
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240 Dif f a=. 1 8 ! DIFF COEF FOR HCL, CM* 2 /S f EST FROM ATER VALUE AND 

DEPENDANCE ON MOLECULAR WT. 

250 Mhcl =36. 47 

260 Mh2o=1 8. 01 528 ! MOLECULAR WEIGHTS 

270 LINPUT "ENTER CASE IDENTIFIER (60 CHRS MAX)", Case? 

280 INPUT "ENTER SURFACE TEMPERATURE DELTA IN CELSIUS = ",Tdelata 

290 INPUT "ENTER DEW POINT DELTA IN CELSIUS = ",Ddelta 

300 INPUT "ENTER WIND DELTA SPEED AT 2 METERS ELEVATION, (m/s) = ",Udelta 

301 Ts=26. 33 

302 Tdew=22. 59 

303 U=. 52 

304 REM INITIALIZE WITH 41 D SPECIFIC DATA 

310 INPUT "ENTER INITIAL ACID CONCENTRATION, WT. PERCENT = ",C 

320 INPUT "ENTER CALCULATION END TIME (minutes) = ",Tf 

330 INPUT "ENTER CALCULATION TIME STEP (minutes) = ",Dtime 

340 INPUT "ENTER INTERVAL FOR EXTRA DATA PRINTOUT (minutes) = ",Dtout 

350 INPUT "ENTER INITIAL LIQUID VOL PER SQ METER =100 CM~3? ",V 

360 INPUT "ENTER MEAN RADIUS OF INITIAL DROPS = 0.1cm? ",Amean 

370 INPUT "ENTER AREA FACTOR FOR EVAPO REDUCTION (4?)",Ra 

380 INPUT "ENTER NAME FOR OUTPUT DATA FILE",Dfile$ 

381 REM FOR Rpt=1 TO 5 

383 REM pt=1 THEN Df ile$="A51 uul " 

384 REM IF Rpt=2 THEN U=2 

385 REM IF Rpt=2 THEN Df ile$="A51 uu2 " 

386 REM IF Rpt=3 THEN U=4 

390 Einf =0 

400 A=0 

410 Cumlgas=0 
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420 Cumlw=0 


430 Sumarea=1 
440 Sigma=.08 

450 Da=. 0075 ! INITIAL BIN SEE FOR RADII IN cm 

460 Tsum=0 ! 

470 GOSUB Prival 

480 REM 
490 REM 

500 REM CALCULATE INITIAL SIZE SPECTRUM 

510 REM 

520 REM 

530 FOR 1=1 TO 40 
540 A( I )=I*Da 

550 Pow(I)=-( (A(I)-Amean) ~2/(2*Sigma~2) ) 

560 Tsum=Tsum+A( I) ~3*EXP(Pow(I) ) 

570 NEXT I 

580 Numo=V/ ( 2 . 0944*Tsum ) I 2/3 PI = 2.0944 

590 Tv=0 

600 FOR 1=1 TO 40 

610 C(I)=C 

620 N ( I ) =Numo*EXP ( Pow ( I ) ) 

630 V(I)=2.0944*N(I)*A(I)~3 

640 Tv=Tv+V ( I ) 

641 PRINTER IS 0 

650 PRINT USING 651 ; I, A( I ) , N ( I ) ,Tv 

651 IMAGE DDD , " A ( I ) =" , D.DDDD, " N ( I ) =" , DDDDDD . DDD , " TOTAL V=",D.DDE 

652 PRINTER IS 16 
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660 NEXT I 


670 REM 

680 REM 

681 REM BEGIN TIME DEPENDANT CALCULATION OF HCL EVAPORATION 

682 FOR T=1 TO Tf STEP Dtime 

683 Ts=26. 33+. 01 05*T+Tdelta 

685 Tdew=22. 59-. 001 9*T+Ddelta 

686 U=.5204+.01 21 *T+Udelta 

711 REM END OF 41 D LAUNCH SPECIFIC DATA 

712 REM 

713 REM 

714 REM 

718 REM CALCULATE VAPOR PRESSURE AND EVAPORATION FOR PURE WATER 
720 REM FLEAGLE AND BUSINGER, P 48 
730 Pwr =9. 4051 -2353/ (Ts+273. 1 6) 

740 Esat=. 750062*1 O'" Pwr ! mm of Hg. 

750 Pwr=9. 4051 -2353/ (Tdew+273. 1 6) 

760 Edew=. 750062*1 O' 1 Pwr 

770 REM PENMAN'S BEST FIT EQUATION ADAPTED TO CURRENT UNITS 

780 REM PENMAN, 1947, PROC. ROY. SOC. A, VOL 193 

790 Lapse=4. 05*60. 0* (1+.526*U)/1 000 

800 Evapo=Lapse* ( Esat-Edew ) ! (g/(m''2 min) see BOOK 8, P9 

81 0 REM 
820 REM 
830 REM 

850 REM 

860 REM 


II-4 



880 


Tsum=0 


890 

Mtg(T)=0 


891 

IF TOO THEN 900 


892 

IF Mtg(T-1 ) >.0001 THEN 900 


893 

Mtg(T)=0 


894 

GOTO 1 460 


900 

II 

o 


910 

Test=0 


920 

FOR J=1 TO 40 


930 

Tsum=Tsum+N ( J ) * A ( J ) 


940 

NEXT J 


941 

IF Tsum=0 THEN 1470 


950 

Kp=Evapo*Sumarea/Tsuin 


960 

FOR 1=1 TO 40 


970 

IF A ( I ) =0 THEN GOTO 1330 


980 

GOSUB Mol 


990 

Lw=2343 !J/g AT 66 deg C 


1000 

Al=25. 624955 


1010 

Bl=. 01 4923 


1020 

Cl=1 .079343 


1030 

La=2343-Al*M( I ) +Bl*EXP (Cl*M ( I ) ) 


1040 

IF M( I ) >8. 06 THEN La=2226 


1050 

REM THESE ARE LATENT HEATS OF VAPORIZATION FOUND BY 

1060 

REM SIMPLE CURVE FIT TO DATA OF A. C. PLEWES 

AT 66 DEG 

1070 

REM ONLY THE RATIO OF LATENT HEATS ENTERS SO 

ONLY ONE 

1080 

REM TEMPERATURE NEED BE CONSIDERED. 


1090 

Khcl=Lw/La*Kp*Dtime 
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1 1 00 Ratio=Dif f w*Mh2o* ( Ew-Eldew )/ (Dif f a*Mhcl* ( Ea-Elinf ) ) 

1110 REM Ratio IS THE WATER TO ACID EVAPORATION RATIO FOR DROP I. 

1120 Delma=Khcl*A ( I )/ ( 1 +Ratio) ! MASS OF ACID LOST FROM Ith DROP. 

1130 Delmw=Khcl*A(I )/( 1 +1 /Ratio) !MASS OF WATER LOST FROM Ith DROP. 

1140 Mass =2 . 0944*A ( I ) ~ 3* ( 1 +. 0049*C ( I ) ) ! INITIAL DROP MASS (GRAMS) 

1150 Massa=Mass*C ( I )/1 00 

1160 Massw=Mass * ( 1 -C ( I ) /1 00 ) 

1170 Newmass-Mass-Delma-Delmw 

1180 IF (1 00*Delma>Massa) OR (1 00*Delmw>Massw) THEN GOTO 1610 
1190 IF Newmass<=0 THEN GOTO 1240 

1200 C ( I ) =(Massa-Delma ) *1 00/Newmass 

1210 A(I)=(Newmass/(2.0944*(1+.0049*C(I) ) ))~ (1/3) 

1220 REM NEW CONCENTRATIN AND RADIUS 

1230 GOTO 1260 

1240 A ( I )=0 

1250 Newmass=0 

1260 V( I )=Newmass/( 1 +.0049*C ( I ) ) 

1270 Mw ( I ) =N ( I ) *De lmw 

1280 Mga ( I )=N( I ) *Delma ! HCL EVAP: g/sq. meter in Dtime 
1 290 Test=Test+Mw( I ) *Lw+Mga ( I ) *La 

1 300 REM 

1310 Mtg ( T ) =Mtg ( T ) +Mga ( I ) ! SUM FOR TOTAL HCL IN TIME INTERVAL 

1320 Tw=Tw+Mw ( I ) ! SUM FOR TOTAL HCL IN TIME INTERVAL 

1330 NEXT I 

1340 Cumlgas=Cumlgas+Mtg(T) 

1350 CumIw=Cumlw+Tw 
1360 GOSUB Env 
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1370 


PRINTER IS 0 


1380 

1390 

1400 

1410 

1415 

1420 

1430 

1440 

1450 

1460 

1470 

1480 

1490 

1500 

1510 

1520 

1 530 

1540 

1550 

1560 

1570 

1571 
1 580 
1590 
1600 
1610 


REM PRINT USING 1 390; T, Mtg ( T) , Tw, Test , Cumlgas ,Cumlw 
IMAGE DDD *DD, 3 { 2X,MD* DDDE ) , 2 ( 2X , DDDD . DDD ) 

IF T MOD 10 THEN GOTO 1420 

PRINT " T MASS HCL MASS WATER TEST 

PRINT USING 1390 ;T,Mtg(T),Tw, Test, Cumlgas, Cumlw 
IF T MOD DtOUt THEN GOTO 1440 
GOSUB Pout 
PRINTER IS 16 
REM 

R EM END OF PRIMARY (TIME) LOOP 

NEXT T 


! CREATE FILE 
I OPEN FILE 


REM 

CREATE Df :T1 4" , 1 4 
ASSIGN Df i le $& " : T1 4" TO #1 
PRINT #1 ?Case$ 

PRINT #1 ; Tf , Dtime , Ts , Tdew , U, C , V, Amean, Ra 
FOR 1=1 TO Tf STEP Dtime 
PRINT #1 ? Mtg ( I ) 

NEXT I 

ASSIGN * TO #1 
REM 

REM NEXT Rpt 
REM 

GOTO 2820 I GO TO END OF PROGRAM 

REM 

REM BEGIN FINE RESOLUTION CALCULATION LOOP 


i CLOSE FILE 


CUM HCL CUM WATER" 


->>>>>>>> 
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1620 REM 


1620 Na=Delma/Massa 

1640 Nw=Delmw/Massw 

1650 IF Nw>Na THEN Na=Nw 

1660 IF Na>1 .5 THEN Na=1.5 

1670 Loop=2+1 34. 0*Na-36. 0*Na*Na 

1680 If Na<0 THEN Loop=121 

1690 Mw ( I )=0 

1700 Mga ( I )=0 

1710 Tivl-Loop 

1720 FOR L=Loop TO 0 STEP -1 

1 730 REM TOP OF SUB LOOP 

1740 IF A(I )=0 THEN GOTO 1300 

1750 GOSUB Mol 

1760 Lw-2343 ! J/g AT 66 deg C 

1770 Al=25. 624955 

1780 Bl=. 014923 

1790 01=1.079343 

1800 La=2343-A1*M(I)+B1*EXP(C1*M(I) ) ) 

1810 IF M( I ) >8. 06 THEN La=2226 


1820 

1830 

1840 

1850 

1860 

1870 

1880 


REM THESE ARE LATENT HEATS OF VAPORIZATION FOUND BY 
REM SIMPLE CURVE FIT TO DATA OF A. C. PLEWES AT 66 DEG C. 

REM ONLY THE RATIO OF LATENT HEATS ENTERS SO ONLY ONE 
REM TEMPERATURE NEED BE CONSIDRED. 

Khcl-Lw/La*Kp*Dtime /Tivl 

Ratio=Dif fw*Mh2o* ( Ew-Eldew )/(Diff a*Mhcl*(Ea-Elinf ) ) 
Delma-Khcl*A( I )/ ( 1 +Ratio) ! MASS OF ACID LOST FROM Ith DROP 
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1890 Delmw=Khcl*A(I)/(1+1 /Ratio) !MASS OF WATER LOST FROM Ith DROP. 

1900 Mass=2. 0944*A(I)~3*(1 +0049*C ( I ) ) ! INITIAL DROP MASS (GRAMS) 

1910 Massa=Mass*C (I )/1 00 

1920 Massw=Mass*( 1-0(1 )/ 100 

1930 Newmass=Mass-Delma-Delmw 

1940 IF Newmass<=0 THEN GOTO 2010 

1 950 C( I ) = (Massa-Delma ) *1 00/Newmass 

1960 IF C(I)<=0 THEN GOTO 2010 

1970 IF C ( I ) >50 THEN GOTO 2010 

1980 A( I )=(Newmass/( 2.0944*(1+.0049*C(I))))~(1/3) 

1990 REM NEW CONCENTRATION AND RADIUS 

2000 GOTO 2030 

2010 A( I )=0 

2020 Newmass=0 

2030 V( I )=Newmass / ( 1 +. 0049*C ( I ) ) 

2040 Mwloop=N( I )*Delmw 

2050 Mgaloop=N ( I ) *Delma ! HCL EVAP: g/sq. meter in Dtiine 

2060 Mw ( I ) =Mga ( I ) +Mwloop 

2070 Mga ( I ) =Mga ( I ) +Mga loop 

2080 Test=Test+Mwloop*Lw+Mgaloop*La 

2090 PRINT USING 21 00; T; I; L; A ( I ) ; Mga ( I ) 

2100 IMAGE "T=" , DDDD.D, " I=",DD," L=" , DDD, " A=" , D.DDD, " 

MGA=",MD.DD 

2110 NEXT L 

2120 REM END OF SUB LOOP 

21 30 GOTO 1 300 
2140 REM 
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2150 REM 


2160 Mol: REM BEGIN SUBROUTINE TO CALCULATE MOLARITY AND VAPOR PRESSURE 

2161 IF C( I ) >40 THEN C(I)=40 

2162 IF C(I)<0 THEN C(I)=. 000001 
2170 M( I ) =1 000*C ( I )/ (Mhcl* ( 1 00-C ( I ) ) ) 

2180 Aw=53.56-8.1475*M(I) + 1 .4352 *M (I)"' 2-. 058161 *M (I)'' 3-1 . 31 7E-5*EXP(M(I ) ) 

2190 Bw=-4. 8693+1 .211 5*M(I)-.21 394*M( I )~ 2+8. 6707E-3*M( I) "3+1 . 981 2E- 

6*EXP (M ( I ) ) 

2200 Cw=-6753.4+363.79*M(I)-66.826*M( I) "2+2.6591 *M ( I )" 3+5. 6202E- 

4*EXP(M( I ) ) 

2210 Aa=30.542-1.3279*M(I)+.03242*M(I)" 2+1 .6501 *LOG(M(I) ) 

2220 Ba=-1 4694+987.61 *M( I )-24.308*M( I)" 2+55.001 *LOG(M(I ) ) 

2230 Ca=9. 1 01 6E5-1 . 58 : E5*M( I )+2045. 9*M ( I )" 2+6331 . 4*LOG(M ( I ) )+ 

23439*M( I ) *LOG(M( I ) ) 

2240 REM ABOVE COEF. FROM DINGLE, NASA CR 2928, JAN 1978. 

2250 REM WATER AND ACID VAPOR PRESSURES ARE OBTAINED AS FOLLOWS: 

2260 Tsk=Ts+273. 16 

2270 Ew=EXP(Aw+Bw*LOG(Tsk)+Cw/Tsk) ! nun OF Hg 

2280 Ea=EXP ( Aa+Ba/Tsk+Ca/Tsk "2 ) ! nun OF Hg 

2290 REM 

2300 RETURN ! ENG OF SUBROUTINE Mol 

2310 REM 

2320 Env: REM SUBROUTINE TO CALCULATE VAPOR PRESSURES NEAR SURFACE 

2330 Sum=0 

2340 Sumarea=0 

2350 Nt=0 

2360 FOR 1=1 TO 40 

2370 Sum=Sum+N ( I ) *A ( I ) 

2380 Sumarea=Sumarea+N ( I ) *Ra*PI*A( I ) *A ( I ) /1 0000 
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2390 Nt=Nt+N ( I ) 

2400 NEXT I 

2410 Abar=Sum/Nt 

2420 Xbar={ 1 0000/Nt) ~(1/3) 

2430 IF Sumarea>1 THEN Sumarea=1 

2440 Ewtemp=Tw/(Lapse*Dtime)+Edew 

2450 Eatemp=Mtg ( T) / ( Lapse*Dtime ) +Einf 

2460 Eldew=Ewtemp*Abar/Xbar 

2470 Elinf =Eatemp*Abar/Xbar 

2480 PRINTER IS 1 6 

2490 IF T MOD 10 THEN GOTO 2510 

2500 PRINT "T " ; "Ewtemp ";"Ew " ; "Eatemp 

" ; "Ea " ; "Abar ";"Xbar" 

2510 PRINT USING 2520;T, Ewtemp, Ew, Eatemp, Ea, Abar, Xbar 

2520 IMAGE DDD, 2X, 4(MD.DE, 2X) ,D.DDDD, 2X,D.DDDD 

2530 PRINTER IS 0 

2540 REM 

2550 RETURN ! END OF SUBROUTINE Env 

2560 REM 

2570 Prival: REM PRINT INITIAL VALUES 
2580 PRINTER IS 0 

2590 PRINT "DATA FILENAME =" ;Df ile$ ,LIN( 2 ) 

2600 PRINT Case$ 

2610 PRINT "SURFACE TEMP (C) =";Ts 

2620 PRINT "DEW POINT TEMP (C) =";Tdew 

2630 PRINT "WIND SPEED (m/s) =";U 

2640 PRINT "INITIAL ACID CONCENTRATION (wt. percent) ="?C 
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2650 PRINT "T final =";Tf; n TIME STEP (min) =";Dtime 

2660 PRINT "INITIAL VOLUME PER SQ METER (cc) =";V 

2670 PRINT "INITIAL DROP RADIUS (cm) =";Amean 

2680 PRINT "AREA FACTOR FOR EVAPO REDUCTION =";Ra,LIN(2) 

2690 PRINT " T MASS HCL MASS WATER TEST CUM HCL CUM 

WATER" 

2700 PRINTER IS 16 

2710 RETURN 
2720 REM 

2730 Pout: Rem PERIOD PRINT OUT SUBROUTINE 

2740 PRINT " I RADIUS VOLUME NUMBER MOLALITY WT PER" 

2750 PRINT " cm cc/drop per sq m " 

2760 FOR 1-1 TO 40 

2770 PRINT USING 2780 ; I, A ( I ) , V ( I ) , N ( I ) ,M ( I ) ,C ( I ) 

2780 IMAGE DD , 2X , D . DDDD , 2X , D . DDDE , 2X , DDDDDD . DDD , 2X , DD . DDD , 2X , DDD . DD 

2790 NEXT I 

2800 PRINT LIN ( 2 ) 

2810 RETURN 

2820 END 
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